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Abstract 

Wc consider the evolution of fronts by mean curvature in the presence of obstacles. 
We construct a weak solution to the flow by means of a variational method, corresponding 
to an implicit time-discretization scheme. Assuming the regularity of the obstacles, in 
the two-dimensional case we show existence and uniqueness of a regular solution before 
the onset of singularities. Finally, we discuss an application of this result to the positive 
mean curvature flow. 
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1 Introduction 

Motivated by several models in physics, biology and material science, there has been a growing 
interest in recent years towards the rigorous analysis of front propagation in heterogeneous 
media, see [271 El US 1211 US] and references therein. In this paper, we analyze the evolution 
by mean curvature of an interface in presence of hard obstacles which can stop the motion. 
Even if this is a prototypical model of energy driven front propagation in a medium with 
obstacles, to our knowledge there are no rigorous results concerning existence, uniqueness 
and regularity of the flow. On the other hand, we mention that the corresponding stationary 
problem, the so-called obstacle problem, has been studied in great detail, see [261 US] and 
references therein. 

To be more precise, given an open set $7 C M", we consider the evolution of a hypersurface 
dE(t), with the constraint E(t) C O for all t > 0, where is an open subset of and M"' \ 
represents the obstacles. The corresponding geometric equation formally reads (we refer to 
Sectional for a precise definition): 

( \ — j '^(^) if X G $7 , , 

"^^^^ ~ \ max(K(x),0) if X € 517 ^ ' 

where v and k denote respectively the normal inward velocity and the mean curvature of 
dE(t) . Notice that the right-hand side of ([T]) is discontinuous on 50, so that the classical 
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viscosity theory |19j does not apply to this case (see however |20l [9] for a possible approach 
in this direction). 

We are particularly interested in existence and uniqueness of smooth (that is C^'^) solutions 
to ([1]). We tackle this problem by means of a variational method first introduced in [5l [23] 
(see also [6] for a simpler description of the same approach), which is based on an implicit 
time-discretization scheme for ([1]). 

After showing the consistency of the scheme with regular solutions (Theorem 14. 8p , we obtain 
a comparison principle and uniqueness of smooth solutions in any dimensions (Corollary 
14. 9|) . Moreover, in the two-dimensional case we are also able to prove local in time existence 
of solutions (Theorem 15. 3p . Notice that in general one cannot expect existence of regular 
solutions for all time, due to the presence of singularities of the flow (even in dimension 2). 
On the other hand, due to the presence of the obstacles, regular solutions do not necessarily 
vanish in finite time and may exist for all times. Eventually, we apply our result to the 
positive curvature flow in two dimensions, obtaining a short time existence and uniqueness 
result (Corollary [63]) for C^'^-regular flows. Indeed, such evolution can be seen as a curvature 
flow where the obstacle is given by the complementary of the initial set. 
We point out that the study of the positive curvature flow in Section [6] is related to some 
biological models which originally motivated our work: in several recent studies of actomyosin 
cable contraction in morphogenesis and tissue repair there is increasing evidence that the 
contractile structure forms only in the positive curvature part of the boundary curve (see 
[U [3] and references therein). Since the contraction of such actomyosin structures can be 
associated with curvature terms (see [22l [H [2]), this leads very naturally to consider the 
positive curvature flow problem. 

Notice that a set evolving according to this law is always nonincreasing with respect to 
inclusion, which is a feature not satisfied by the usual curvature flow. This shows why 
assembling the contractile structure only in the positive curvature portion of the boundary 
(instead of all around) and thus doing positive curvature flow (instead of usual curvature 
flow) is an interesting way to evolve from the biological point of view: it corresponds to 
making our wound (or hole) close in a manner where we never abandon any portion of the 
surface we have already managed to cover since we started closing. 

We also remark that the positive curvature flow is useful in the context of image analysis [28t 
p. 204], and appears naturally in some differential games |23j . 

2 Notation 

Given an open set A C M", a function u G L,^{A) whose distributional gradient Du is a 
Radon measure with finite total variation in A is called a function of bounded variation, and 
the space of such functions will be denoted by BV{A). The total variation of Du on A turns 
out to be 

supjy" u divz dx:z£ C^{A;W), \z{x)\ < 1 Vx e ^| , (2) 

and will be denoted by |Du|(A) or by \Du\. The map u — )• |Dn|(A) is L-'^(j4)-lower semi- 
continuous, and BV{A) is a Banach space when endowed with the norm ||ii|| := \u\ dx + 
\Du\{A). We refer to [7] for a comprehensive treatment of the subject. 
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We say that a set E satisfies the exterior (resp. interior) i?-ball condition, for some R > 0, 
if for any x G dE there exists a ball Br{x'), with x G dBR{x') and Bji{x') f] E = ^ (resp. 
Br{x') C E). Notice that a set E with compact boundary satisfies both the interior and the 
exterior i?-ball condition, for some > 0, if and only if dE is of class C^'^. 



3 The implicit scheme 

Following the celebrated papers [5l [23] , we shall define an implicit time discrete scheme for 
([1]). As a preliminary step, we consider solutions of the Total Variation minimization problem 
with obstacles; the scheme is then defined in Definition 14.21 below. 

Let B C M" be an open set and let f : i? — )• [— oo, oo) be a measurable function, with G 
L2(S). Following [ana [Ml, given > and / G L'^{B), we let Sh^y{f,B) G L'^{B)nBV{B) 
be the unique minimizer of the problem 

min / \Du\ + \ [ {u - ff dx. (3) 

u>v 2h Jb 

We have the following comparison result (see \15\ Lemma 2.1]). 

Proposition 3.1. The operator Sh, {-,B) is monotone, in the sense that ui = Sh,vi{fi, B) > 
U2 = Sh,v2{f2, B) whenever fi > I2 and vi > V2 a.e. 

Proof. The idea is simply to compare the sum of the energies of ui and n2, with the sum of 
the energy of ui A U2 (which is admissible in the problem defining U2) and of ui V U2 (which 
is admissible in the problem defining ui). The conclusion follows from the uniqueness of the 
solution to ([3|). □ 

Proposition 3.2. Assume f,v+ G L°°{B): then u = ShM^B) G L°°{B) and 

\\Sh,v{f^^)\\L°-{B) < max (||/||loo(b), \\v^\\l^(b)) ■ 

Proof. Again, the proof is trivial. It is enough check that the energy of um = {u\/ —M) /\M is 
less than the energy of while um is admissible as soon as M > max Ib^ IIl°°(b)) • 

□ 

Theorem 3.3. Let v : R"' ^ [-co, +00) be a measurable function with G L^^{W^), 
f G L[^^(M"), and h> 0. There exists a unique function u G L^^{W^) Ci BV'iocO^"') , which we 
shall denote by Sh,v{f), such that for all R > and p G (n, +00) there holds 

M ™oo ~ ^h,v{f, BM)\\LP{Bn) = 0- 

This function is characterized by the fact that u > v a.e., and for any R and any ^ G i?y(M") 
with support in Br and u + (p > v a.e., 
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Proof. We shall show a bit more: for any M > 0, let us denote by um an arbitrary local 
minimizer of (l3|), in the sense that 

for any G BV{Bm) with compact support. We will show that {um)m>2R is a Cauchy 
sequence in LP{Bji), provided p > n. 

To start, let us consider ip : M ^ M_|_ a smooth, nondecreasing and bounded function with 
< ^p{s) < Cs+ for any s. Let M' > M > 0, and let ip G C;?°(Sm; which we extend by 
zero to Bm'- We denote u = um, u' = um'- Let t > 0: observe that 

u' {x) + til){u{x) — u' {x))ip{x) > u'{x) > v{x) 

u{x) — t'ip{u{x) — u'{x))ip{x) > u{x) — tC sup if {u{x) — u' {x))~^ 

> u{x) — {u{x) — u'{x))~^ 
= mm{u{x) , u' {x)} 

> v{x) 

for almost every x G M", as soon as t < (C sup^?)""*^. 
Hence, we deduce from (jl]) that for t small enough. 



\D{u — ttlj{u — u')(p)\ + — I \u — tilj{u — u')(p — f\'^ dx 



> I \Du\ + ^ [ \u-f\^dx 
JBm -JBm 



and 



\D{u' + til;{u-u')ip)\ + ^ / \u' + til^iu - u')ip - fl"^ dx 
Bm ^'^JBm 



> [ \Du'\ + ^ f \u'-f\^dx., 
JBm ^'^ JBm 

which we sum to obtain 

. {u — u')ip{u — u')ip dx < — {il){u — u')ip)'^ dx 
h Jbm Jbm 

+ / \Du-t'il)'{u-u ){Du- Du')ip-ti;{u-u')Vip\ 



Ba 

+ \Du' + tip'iu - u'){Du - Du')ip + t^{u - u')V(p\ - \Du\ - \Du'\ . 

For p < t\\ip\\cxi\\ip'\\oo ^ 1 and t small enough, the integrand in the right-hand side has the 
form 

\p - p{p - v') - + \v' + p{v - p') + iq\ - \p\ - \q\ 

< 2t\q\ + {l-p)\p\+p\p'\ + il-p)\p'\ + p\p\-\p\-\q\ = 2t\q\ 
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and we obtain 



t 
h 



{u — u')il){u — u')ip dx < — {il){u — u')Lpf' dx + 2t I iIj{u — u')\Vip\dx . 

h .In., Jb,, 



iBm JBm 
Dividing by t and letting t — >• 0, we deduce 



{u — u)'il){u — u')ip dx < 2h ip{u — u')\V(p\dx . 



(5) 



Consider now, for p > 2, the function ip{s) = (■s"'")^"^: we want to show that (0) still 
holds. We approximate ip with ipk{s) = fetanh('0(s)/A;), for k>\. The functions ipk satisfy 
the assumptions which allowed us to establish ([5D, so that it holds with ip replaced with 
ipk- Moreover, \\m.k^oo i^kiu — u') = sup^>;^^fc(ti — u') = '4'{u — u'), and in the same way 
supfc>i(u — u')il)k{u — u') = {u — u')'ip{u — u'). Hence, the monotone convergence theorem 
shows that ([5]) also holds, in the limit, for '(/', as claimed. 

We can take ip{x) = 'po{\x\/MY, for some G C~([0,1);M+) which is 1 on [0,1/2]. It 
follows from ([5]) and Holder's inequality that 



[{u-u')+M\A/M)Y' dx 



< 2h I [{u - uYM\^\/M)Y-' I- \^',{\x\/M)\ dx 
'Bm 



< 2h 



Bk 



[{u-uYM\x\/M)y 



i-i 

p 



Ba 



P \P 

M 



WM/M)\ 



Hence: 



(u-u )+^o ( -j^ 



< 



2hpul/^ 

]\^l-n/p I 



\fo\ 



Lp(Bm) 

with ijjn the volume of the unit ball. Exchanging the roles of u and u' in the previous proof. 



we find that 



I MM - '^M'\\lp{Bm/2) 



< 



2hpuin^ 



(6) 



As in particular um (or um') could, in this calculation, have been chosen to be the minimizer 
Sh.vifi Bm), which is bounded by Proposition 13.21 we obtain that um' & (as well 

as Um)- Hence, choosing i? > 0, we see that {um)m>2R defines a Cauchy sequence in I/'{Bji), 
provided p > n. It follows that it converges to some limit u S LP{Bji). As R is arbitrary, we 
build in this way a function u which clearly satisfies the thesis of the theorem. □ 



Corollary 3.4. Assume f> f , v > v' , h > 0, then Sh,v{f) > Sh,v'{f')- 
Proof. It follows from Proposition 13. II and the definition of Sh^vif)- 



□ 



Corollary 3.5. If f,v are uniformly continuous on R", with a modulus of continuity u}{-), 
then Sfi^j;{f) is also uniformly continuous with the same modulus of continuity. 
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Proof. It follows from the previous corollary. For z € M", let := v{x — z)—uj{\z\) < v{x) 

and fix) := f{x - z) - uj{\z\) < f{x). Then, Sh,,'{f') = - z) - u;{\z\) < S^Af), 

which shows the corollary. □ 

Observe that, if f,v are uniformly continuous, then Sh^vifiB) satisfies the elliptic equation 

u - f 

— div z H ; — =0 on {x G i? : u(x) > v(x)}. (7) 

h 

where the vector field z satisfies |2;| = 1 and z = Du/\Du\ whenever \Du\ ^ 0. 

Proposition 3.6. Assume that f{x) —)• oo as |rE| — >• oo, and let s G M. Then the set 
{Sh,vif) < s} is the minimal solution of the problem 

min P(E) + / dx. (8) 

Ec{v<s} Je h 

Similarly, the set {Sh^vif) ^ s} is the maximal solution of 

min P(E) + [ dx. (9) 

Ec{v<s} Je h 

Proof. Let M > and consider the set E^^ = {Sh,v{f, Bm) < s}. Reasoning as in [11] (see 
also Sec. 2.2.2]) one can show that Ef_j is the minimal solution of 



min P{E,Bm)+ [ ^^dx. 
tMn{v<s} Je h 



EcBMr\{v<s} 

Since / is coercive, the sets Ef^.j- do not depend on M for M big enough, and coincide with 
the set {Sh,vif) < s}, so that the result follows letting M — > +oo. 

The second assertion regarding the set {Sh,vif) < s} can be proved analogously. □ 

4 Mean curvature flow with obstacles 

Let us give a precise definition of the flow ([1]) . Given a set E C M" we denote by 

dEix) := dist(x, E) - dist(x, W\E) x G M" 
the signed distance function from E, which is negative inside E and positive outside. 
Definition 4.1. Given a family of sets E{t), t G [0, T], we set 

d{x,t) := dE(t){x). 

We say that E{t) is a C"^'^ supersolution of ([1]) if there exists a bounded open set U C M" 
such that E{t) C and dE{t) C U for all t G [0, T], 

d G Lip([/ X [0,T1) 

(10) 
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and 

dd 

— > Ad + 0{d) a.e. in [/ x [0, T]. (11) 
We say that E{t) is a C^'^ suhsolution of ([T]) if pT|) is replaced by 
dd 

— <Ad + 0{d) a.e. in {U x[0,T])n{d> dn}, (12) 

and we say that E[t) is a C^'^ solution of ([1]) if it is both a supersolution and a suhsolution. 

We now fix an open set C M" (representing the complement of the obstacle) and a compact 
set £^ C fi. The case when E'^ is compact can be treated with minor modifications. 
Since E is compact, without loss of generality we can assume that O is bounded. Indeed, as 
it will be clear from the sequel, replacing il. with n B^j will not affect our construction, 
provided Bj^j D E. 

Definition 4.2. Let h > and set 

ThE:={Sh,d^idE)<0}. (13) 

Given t > 0, we let 

Eh{t) := rf^'i? 
he the discretized evolution of E defined by the scheme Th. 

Notice that T^E is an open subset of 17 and, by Proposition 13. 6t T^E is the minimal solution 
of the geometric problem 

minP(F) + i / dsdx (14) 
Fcn ^ ' h Jp 

or equivalently 

minP(F) + ^ / \dE\ dx. 

FCQ hJp^E 

When n = M" this corresponds to the implicit scheme introduced in [5l [23] for the mean 
curvature flow. Here, from ^ it also follows that T^E satisfies 

K + ^ = on dThE \ dVL. (15) 

Remark 4.3. Observe that from Proposition 13. II it follows 

EiCE2^ ThEi C ThE2. 

Moreover, by Corollary 13.41 we have Sh,dci{dE) > Sh^^oo{dE) which implies T^E C T^E := 
{Sh,~oo{dE) < 0}. Notice that T^E is the scheme introduced in [SI [23] for the (unconstrained) 
mean curvature flow. 

From the general regularity theory for minimizers of the perimeter with a smooth obstacle 
, [T2] we have the following result. 
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Proposition 4.4. Let be of class C^'^, E Q Q and h > 0. Then there exists a closed 
set S C dThE n Q such that = for all s > n - 8, dT^E \T. is of class C^'S and 

{dThE nn)\^ is C^'" for any a < I. 

Proposition 4.5. Let dVL he of class C^'"^ . Then there exists C{Q) > such that 

ThE = {Sh,-oc{dE + Chxnc) < 0} 
for all C > C{Vt). In particular T^E is a minimizer of the prescribed curvature problem 

inmP{F) + C\F\n\ + ^ [ dsdz. (16) 
F h Jp 

Proof. We recall that Sh,-ooidE + Chx^c) is the limit, as M — )• oo, of the minimizer um of 
the variational problem 

min / l-Dnl H / (u — dp — Chxoc)'^dx (17) 

ueBViB,j)JB,, ^hJB,, 

From Proposition 13.61 it follows that T^E is the minimal solution to (|14p . while 

F = {Sh,^oo{dE + Chxnc) 

is the minimal solution to (jl6p . li F d Vt, then |F \ J7| = and both F and T^E solve the 
same problem, and they must therefore coincide. 

In order to show that F C $7, it is enough to find a positive constant C such that for all 
X ^Q,, Um > C > for M large enough. 

By assumption, Q. satisfies an exterior i2-ball condition, for some i? > 0, that is, for any 
a; ^ fi, there is a ball B^ix') with x G Bfi{x') and Bb{x') R $7 = 0. If M is large enough, we 
also have Bpi{x') G i?A//2- Since E d dp + hCxn^ > hCxsn^x')-, so that um is larger than 
the minimizer u' of 



/ l^^l + ^ / ~ hCxBR{x')f dx 



mm 

u&BV{Bm) Jbm ' ' ''^'^ Jb 



If C > n/R, then it is well known that for M large enough, u' > (C — n/R)h a.e. in 
XBr(x') ES]. The thesis then follows. □ 

4.1 Existence of weak solutions 

As a consequence of Proposition [531 when dO. is of class C^'^ the scheme enters the framework 
considered in |17] . In that case, we can also show existence of weak solutions in the sense 
of O [25. We observe that the results in [6l p. 226] still apply and we can deduce the 
(approximate) l/(n + l)-Holder-continuity in time of the discrete flow starting from an initial 
set Eq. As a consequence, following [U Th. 3.3], we can pass to the limit, up to a subsequence, 
and deduce the existence of a flow Fit), which is Holder-continuous in time in L^{Q). 

Theorem 4.6 (Existence of Holder-continuous weak solutions). Let dO. be of class C^'^, 
let E d be a compact set of finite perimeter and such that \dE\ = 0. Let Eh{t) be the 
discretized evolutions starting from E, defined in Definition \4-S\ Then there exist a constant 
C = C{n,E,Vt) > 0, a sequence hi ^ and a map E{t) "Pi^) such that 
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• E{0) = E; 

• E{t) is a compact set of finite perimeter for all t > 0; 

• linij \Ehi{t)AE{t)\ = for all t > 0; 

• \E{t)AE{s)\ < C\s - for all s,t>0, with |s -t| < 1. 
4.2 Consistency of the scheme 

The main result of this section (Theorem l4.8p is showing that the implicit scheme is consistent 
with regular evolutions, according to the following definition. 

Definition 4.7. The scheme is consistent if and only if 

1. If E{-) is a supersolution (see Def. \4-l^ in ctn interval [ti,t2], then for any t € [^1,^2]; 
any Hausdorff limit of TJ^E{ti), n — )• 00, /i — > 0, nh ^ t — ti, contains E{t). 

2. If E{-) is a suhsolution, this inclusion is reversed. 
Theorem 4.8. The scheme T^ is consistent. 

Proof. The proof consists in building, arbitrarily close to dE{t), strict super and subsolutions 
of class C^, of the curvature flow with forcing term Cxn<=, for C large enough. Then, the 
consistency result in [T71 Th. 3.3] applies. 

Step 1. Let E he a subsolution on [ti,t2] in the sense of Definition 14.11 let U C M" be the 
neighborhood associated to dE{t) (given by Definition 14. Without loss of generality we 
can assume ti = 0. 

Observe that there exists p > such that < p} C f7 for all t G [0,t2]> and the sets 

9r2, d{d{-,t) < s}, \s\ < p, satisfy the interior and exterior p-ball condition for all times (in 
particular dE{t) satisfies the condition with radius 2p). 
Let Cp > {n — 1) / p'^ , and for e > small, let 

ds{x, t) = d{x, t) — e — AcpEt t G [0, t2]. 

Observe that for e small enough, {|de(-,i)| < p/2} C < p} for all t. The constant 

Cp is precisely chosen so that in this set, the curvature of two level surfaces {d{-,t) = s} and 
{d{-,t) = s'} at points along the same normal vector V(i(-,t) differ by at most Cp\s — s'\. 
We have, for a.e. t G (0,t2) and x G {\d{-,t)\ < p} CU, 

dde / N 9d ,^ / N N , 
-^(a;,*) = -Q^{^dEix,t)[x),t) -4:Cpe, 

thus: 

• If IlQE(x,t){x) G ^, then (by Definition 14. ip 

dd 

-Q^i^^t) < Ade{x,t) -Acpe + Cp\d\ < A4(x, t) + Cpj4| + Cp(-4e + e(l + 4cpt)) 

so that iit<t = min(t2, l/(2cp)) and {d^l < e/2, 

dd s 

^(X,t) < Ade{x,t) - Cp-. (18) 
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• While if IlQ£;(x,t)ix) S on, then d = and almost surely dd/dt = 0, so that dde/dt = 
—4:Cp£. On the other hand, there is a constant C large enough (of order 1/p, and 
admissible for Proposition 14. 5p such that lAd^l < C a.e. in < p}, and we 

deduce 

dd 

-AcpE = -Q^ix,t) < Ads{x,t) + C -4cp£. (19) 

Moreover, if d^ > — e/2, we have that dn = d > AcpEt + e/2. 

Consider a function gg, which is C in {do, > e/2}, in il, and smoothly decreasing from C to 
as dn decreases from e/2 to 0: we deduce from (jlSp and (jl9p that 

(9(ie . , e 

a.e. in {(x, : \d^{x, t)\ < e/2 , t G (0, f)}. We have built a strict subflow, as close as we want 
from dE{t), for i € [0, t]. The fact that t could be less than t2 is not an issue, as we will see 
in the end of the next step. On the other hand, the consistency result in |17j requires that 
d is at least in space, which is not the case here (and the proof does not extend to C^'^ 
regularity). For this, we need an additional smoothing of the surface, which we perform in a 
second step. 

Step 2. Now consider a spatial mollifier (f-qix) = r]~'^{p{x/rj), with rj « e. For all time let 
de = * de, which is still Lipschitz in t and now, smooth in x. If r] is small enough, and 
since gg is continuous, we have 

ddg ^ ,„ e 
— < A4' + ..-c,- 

for a.e. x,t with |(i£(a;,t)| < e/2 — ij. We can rewrite this equation as a curvature motion 
equation with some error term, as follows: 

Now, we have that 

1 > jVd^l > I-CT? (21) 

almost everywhere, for some constant c > 0, of order 1/p. Hence, if r] is small enough, we 
have 

ge{l-\Vd'^\)<Cpe/16. (22) 

We claim that the following estimates holds: there exists a constant c > (of order 1/p^) 
such that 

p^d^Vd^l < C7]. (23) 
This will be shown later on (see Step 3). Using ([2T]l and ([23]) . we find that 

— WW — - 
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if T] is small enough. Thus (|20p becomes, using ()22p . 



Since |-D^(ie| < 1/p for a.e. t and x with |(f£(x,t)| < e/2, this is also true for \D'^cQ\ (for 
< e/2 — 77), and using ([2T]l we can easily deduce that the boundaries of the level 
sets E^{t) = {d£{-,t) < 0} have an interior and an exterior ball condition with radius p/2. 
Together with (|24|) . and using < Cxw, we find that Ei;{t), < t < t, is a strict subflow 
for the motion with normal speed V = —k — Cxfi'^, and [T71 Th. 3.3] holds. We deduce that 
there exists ho > such that if /i < /iq, Th{Eir{t)) C Efr{t + h) for any t £ [0,i — h], where 
is the evolution scheme defined by 

Tf,E = {Sh,-oo{dE + Chxn^) < O} 

for any bounded set E. (It corresponds to the time-discretization of the mean curvature flow 
with discontinuous forcing term — Cxn=-) Recall that if E C il, Proposition 14.51 shows that 
TfiE = TfiE C 0. In particular, for the subflow E{-) considered here, he have TJ^{E{0)) = 
Tf^E{0), for all n and h > 0. By induction, it follows that as long as nh < i, 

TJiE{Q) = riE{0) C E,{nh), 

hence Tj;^''^E{0) is in a 3e- neighborhood of E{t). Since t only depends on p > (the 
regularity of the subflow E{-)), we can split [0,t2] into a finite number of intervals of size at 
most t and reproduce this construction on each interval, making sure that the e parameter 
of each interval is less than one third of the e of the next interval. 

We deduce that for any (5 > 0, if /i > is small enough, then TJ^E{0) C {dE[nh) ^ foi' 
< nh < t2- This shows the consistency of with subflows, assuming (|23|) holds. 

Step 3: Proof of estimate ([23|1 . Recall that since de is a distance function, iVd^l = 1 almost 
everywhere. Now, let us compute, for 77 > small and a;,y € {d{-,t) < e/2 — r/}: 

|Vff^(x,t)|2 - |Vd^(y,t)|2 = (Vd^(x,t) - Vd^(y,t)) • (Vd^(x, t) + Vd^(y, t)) 




{Vde{x - z,t)- Vdeiy - z,t))-{Vde{x - z',t) + V4(y - z' ,t)) ipr,{z)ipr,iz') dz dz' . 



(25) 

As \D'^de\ < 1/p, Vds{-,t) is 1/p-Lipschitz, using |V4(a; - z,t)\'^ - |V4(y - = it 

follows 

{Vdeix -Z,t)- V4(y - Z, t)) ■ {Vdeix - z' , t) + V de{y - z' , t)) 

2 2 

< \Vde{x — z,t) — Vdeiy — z,t)\-\z — z'\ < —^\x — y\\z — z\ 

P 

and it follows from (|25p that 

|V4^(x,t)|2-|Vd?(y,t)|2 < \\x-y\^. 
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We deduce (letting y — x) that 

2\D^d'J{x,t)Vd'J{x,t)\ < ^rj, 

which is estimate (j23|) . 

Step 4- Consistency with superflows: the proof is almost identical (reversing the signs and 
inequalities), but simpler for superflows. Indeed, all the sets we now consider stay in Q and 
we do not need to take into account the constraint or the forcing term Cxw- 

□ 

We can define a generalized flow as limit of the scheme as /i — > 0. Given an initial set 
E cn, for all t > we let 

Eh{t)=TfJ''^E and E^ = [j Eh{t) x {t} dW' ^ (26) 

t>Q 

Then there exists a sequence {hk)k>i such that both E^^ and M" x [0, +cx)) \ E^^ = '^E^^ 
converge in the Hausdorff distance (locally in time) to E* and ^E^ respectively. 
From Corollary 13.41 and Theorem 14.81 we obtain a comparison and uniqueness result for 
solutions of dl]). 

Corollary 4.9. Let Ei(t) and E2{t) be respectively a sub- and a supersolution of ([TJ for 
t G [0,T], in the sense of Definition \4- 1[ Then, if Ei{0) C £^2(0), it follows that Ei(t) C E2{t) 
for all t S [0,r]. In particular, if dE is compact and of class C^'^, there exists at most one 
solution E{t) starting from E. Moreover, by Remark \4-3[ E[t) is contained in the solution 
to the (unconstrained) mean curvature flow starting from E. 

5 Short time existence and uniqueness in dimension two 

In this section we assume n = 2 and dO. of class C^'^. In the bidimensional case, the mean 
curvature is the same as the total curvature of the boundary dE. Hence, any estimate on the 
mean curvature yields a global estimate on the regularity of E. This will be the key of our 
construction, for showing the existence of regular (C^'^) solutions to the mean curvature flow 
with obstacles. In higher dimension, this is not true anymore, and showing the existence of 
such solutions remains an open problem. 
The following result follows as in |lll Lemma 7]. 

Lemma 5.1. Let h > and let E Q, with dE of class C^'^. Let 6e be the maximum 6 > 
such that both dE and di} satisfy the 5-ball condition, and let u = Sh,dn{dE)- Then, for all 
5' S (0, 6e) we have 

\u-dE\<^^ m{\dE\<5'} (27) 
Oe - 

for all h<{6E- 

Lemma 5.2. Let ECU with dE of class C^'^. Then, there exists 5 > and T > such 
that 

dEh{t) satisfies the 5-ball condition for all t G [0,T]. (28) 
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Proof. Let 5e he as in Lemma l5.H and let K = 2/5e- By Lemma l5. 11 applied with 6' = Kh, 
we get 

h h / K, ^K"^ 



inmE. BE) < it^^ + f/^ C--^/ 

for all h < ho := 51^/12, where the constant C > is independent of E. Recalling (jlSp and 
Proposition 14.41 we get 

„ „ 1 / iT, p,K'^ ,^ 



which implies 



5t^e > minf^ ,6e - dn{dThE,dE)) (29) 

\\\l^\\L°°idThE) J 

( h ( K ^K"^ \ f K ^K'^ 
> • min 1 - 1 + —h + C-^h^ , 1 + -^h + C-^h' 



y 6l \ 6e Si J ' \ 6e 6\ 
for all /i < /iQ. By iterating ^ we obtain □ 
We now prove a short time existence and uniqueness result for solutions to ([T]). 

Theorem 5.3. Let dO. be of class C^'^ and let E Q Q. with dE of class C^'^. Then there 
exists T > such that dTJ admits a unique C^'i solution E(t) on [0, T] with E{0) = E. 



Proof. Let E^ he as in (p6|) and let 





't' 


a- 


h 



^Eh{t+h) 



By Lemmas 15.11 and 15.21 there exist an open set U C M" and T > such that dE]^[t) C 
U for all t G [0, T] and \V'^dh\ € -L°°(C/ x [0,r]); moreover, recalling ([27D we also have 
d/j € Lip([/ X [0, T]). By the Arzela-Ascoli Theorem the functions d^ converge uniformly in 
U X [0, r], up to a subsequence as ^ 0, to a limit function d G Lip(?7 x [0, T]) such that 
iV^d] G L°°(;7 X [0,r]) and \Vd\ = 1 in C/ x [0,T]. Letting E{t) = {x : d{x,t) < 0}, for all 
t G [0,r] we then have E{0) = E, E{t) C 9. and dE{t) is of class C^'^ 

It remains to show that (jlip and ()12p hold in [/ x [0,T]. From Theorem 14.81 it follows that, 
given a supersolution E{t) on [^1,^2] C [0,T] with E{ti) C E{ti), we have -E(i) C for all 
t G [ii,i2]i and the same holds with reversed inclusions if E{t) is a subsolution. This implies 
that 

dd 

— = Ad a.e. in {U x [0, T]) n {d > d^} n {d = 0}, 

which proves (|12p . Observe that, by parabolic regularity, dE[t) n is an analytic curve and 
the equality holds everywhere. 
As we have 

dd 

— = a.e. in (i7 X [0,r]) n{d = df^}, 
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the proof of (jlip amounts to show 



Ad < a.e. in {U x [0, T]) n {d = dn}. (30) 
Assume by contradiction that there exist (x, t) E ([/ x {0,T)) {d = d^} such that 

dd 

— {x,t) = 0<Ad{x,t) = Adnix). (31) 
ot 

Without loss of generahty we can assume d{x,t) = dn{x) = 0, and d^ is twice differentiable 
(in the classical sense) at x. 

Let us take an open set J7 D with (compact) boundary of class C°° and such that 

X e dfl and ^d^{x) > AdQ{x) > 0. 

We let r2(t), for t G [0, r] and r > 0, be the evolution by curvature of [5j, and observe that 
E{t) = i}{t — t), t (z [t,t + T], is a subsolution in the sense of Definition 14.11 In particular, by 
Theorem US] 

E{t) C E{t) for all t G [t, t + r], 
but this implies, letting d{x,t) = d^^^^{x) and recalling ([HT]) . 



-(x t) > -! 

leading to a contradiction. This proves ()30p and thus ()lip . 

Finally, the uniqueness of E{t) follows from Corollary 14.91 □ 

Remark 5.4. Notice that in Theorem 15.31 it is enough to assume that il. satisfies the exterior 
i?-ball condition for some i? > 0, which is a weaker assumption than requiring d^l to be of 
class C^'^. Indeed, we can approximate Q with the sets 

^p-= U ^p(^) p>^- 

Bp{x)cn 

Notice that C [7 and SJlp is of class C^'^, for all p > 0. If we take p small enough so that 
E 0,p then, by Theorem 15.31 applied with replaced by 0,p, we obtain a solution Ep{t) on 
[0,Tp]. However, Ep{t) is also a solution of the original problem, with constraint 17 instead 
of Op, since Jlp is a subsolution to ([1]) in the sense of Definition 14.11 



6 Positive mean curvature flow 



In this section we consider the geometric equation 

f = max(K, 0). (32) 

Notice that, by passing to the complementary set, (|32p includes the evolution by negative 
mean curvature v = min(K, 0). 
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Definition 6.1. Given a family of sets E[t), t € [0, T], we set 

d{x,t) := dE{t){x)- 

We say that E{t) is a C^'^ solution of ([32]) if there exists a bounded open set U C M" such 
that dE{t) C U for all t G [0, T], 

d G Up{U X [0, T]) I V^d] GL'^iU X [0, T]) 

and 

dd 

— = max (Ad, 0) + 0{d) a.e. in [/ x [0, T]. (33) 

Lemma 6.2. Lei -Ei(t) and E2{t), with t £ [ti,t2], two C^'^ solutions of ([32]) m t/ie sense 
of Definition \6.1[ Then, if Ei[ti) C E2{ti), it follows that Ei{t) C i?2(i) for all t S [ti,t2]- -^^ 
particular, if dE is compact and of class C^'^ , there exists at most one solution E(t) starting 
from E. 

Proof. Notice that it is enough to prove the thesis with t2 = + for some r > 0, since the 
general claim then follows by iteration. Fix e > 0, let de(x,t) := dE2(t){x) + £ + Cet and let 

Ee{t) := {x : de{x,t) < 0} t G + t], 

where the positive constants C, r will be determined later. Notice that dE^ is compact and 
of class C^'^ for all e small enough, and i?e(t) E2{t) as e — )■ 0. A direct computation gives 



dt 
where 



dd 

' > max {Ade + e{C- CK^t - K^) , O) + 0(4) a.e. inU x [ti,ti + r], (34) 



K= sup \\AdE2{t)\\L^{dE2(t))- 

X&[tl,t2] 



If we choose C = 2K'^ and r < 1/C, (fM|) implies that -Ee(t) is a supersolution of ([32]) . 
Letting -De(t) := dist(9£'i(t), 9E'e(i)), we have that -D^ is Lipschitz continuous, -0^(0) > e 
and ^^(t) > for a.e. t € [ti, ti + r]. As a consequence, Eir{t) C £^i(t) for all t € [ti,ti + t], 
and the thesis follows by letting e — ?• 0. □ 

Remark 6.3. Notice that the viscosity theory [TO] applies to equation (j32|) . since the function 
K —7- max(K,0) is continuous. Then, Lemma 16.21 implies that, if the initial set has compact 
boundary of class C^'^, the corresponding viscosity solution does not create fattening, i.e. is 
unique, before the onset of singularities. Corollary 16.51 below will establish the existence of 
such C^'^ solutions. 

Given E C M" and h > 0, we set E^ = E^ = E and, by iteration, 

for ah n G N. We also let Eh{t) := eI'^"^ and Eh{t) := eI^^\ Notice that Eh{t) is the 
discretized evolution corresponding to the mean curvature flow with obstacle VL = E (see 
Definition 14. 2p . while Eh{t) is an implicit scheme for (j32p . 



(35) 
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Proposition 6.4. Let h > and let E C he a set with compact boundary. Then 

Eh{t) = Eh{t) for all t>0. 

In particular 

Eh{t2) ^ Eh{ti) for all ti < is- (36) 

Proof. We have to show that E'^ = E'^ for all n G N. By the definition we have Ej^ = E^ =: F. 
If we also show that Ef^ = £^^, then the thesis follows by iteration. As di? > d^;, by Proposition 
13.11 we have that Sh^d^idF) > Sh,dEi'^E), so that 

El = {Sh,dEidF) < 0} C {Sh,d^{dE) < 0} = F. (37) 

By Proposition 13.61 we know that E^ is the minimal solution of 

min P(X) + \ I dpdx. 
XCE hJx 

Recalling ()37p it then follows that E^ is also the minimal solution of 

min P(X) + i / dpdx 
XCF hJx 

and hence coincides with Ef^, again by Proposition 13.61 □ 

Proposition 16.41 implies that the evolution ()32p . with initial set E, can be seen as a particular 
case of ([1]) with ft = E. As a consequence, from Theorem 15.31 we get a short time existence 
result for regular solutions to ()32p . 

Corollary 6.5. Let E C M."^ with compact boundary of class C^'^. Then there exists T > 
such that (|32p admits a unique solution E{t) on [0, T] with E{0) = E and dE{t) a compact 
set of class C^'^ for all t G [0, T]. Moreover 

E{t2) C E{ti) for all ti < tg- (38) 

Proof. Thanks to Theorem 15.31 there exist T > and a unique solution E{t) of ([T]) on 
[0,r], with E{0) = E = n and dE{t) of class C^-i. By Proposition [631 for ah t G [0,r), 
E{t) is the solution of ([1]) on [t,T] with obstacle = E[t). In particular, letting as above 
d{x,t) = (i£;(j)(x) and recalling (fTT]l . this implies 

dd 

— = max{Ad + O{d),0) a.e. in [/ x [0, T] , 
that is, E{t) is the solution of (j32p in the sense of Definition 16.11 

The uniqueness of E{t) follows from Lemma 16.21 and ()38p follows from (|36p. □ 
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